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We examine the possibility of achieving quintessential inflation, where the same field serves as both 
infiaton and quintessence, in the context of a five-dimensional braneworld. Braneworld cosmology 
provides an appropriate environment as it permits inflation with much steeper potentials than the 
conventional scenario, which is favourable to a late-time quintessence. We explore a wide space 
of models, together with contemporary observational data, to determine in which contexts such a 
picture is possible. We find that such a scenario, although attractive, is in fact impossible to achieve 
for the potentials studied due to the restrictiveness of current data. 

PACS numbers: 98.80.-k 
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I. INTRODUCTION 

In the past decade there has been a lot of interest in 
the study of extra dimensions, in particular the scenario 
where our apparent Universe lives in a three-dimensional 
hypersurface, a 3-brane, embedded in a space with more 
dimensions (for a review see Ref. pQ). This is the so-called 
braneworld picture, which brings new possibilities and 
perspectives to cosmology. In this work we will consider 
the simplest case of a braneworld scenario, where there 
is only one extra dimension, and analyze the impact of 
this modification in the attempt of getting quintessential 
inflation [2]. 

Our framework is of Kaluza-Klein type, but where, 
instead of compactifying the extra dimension, one as- 
sumes that matter fields are confined to the brane and 
that gravity is allowed to act in the five-dimensional em- 
bedding space, generally called the bulk. We adopt the 
simplest such set-up, known as Randall-Sundrum Type 
II [3] . By solving the Einstein equations in the vicinity of 
the brane, one finds a modified Fricdmann equation that 
determines the cosmological evolution caused by fields 
confined to our 3-brane, such as a scalar field [4]. In 
this case, the inclusion of an extra dimension in the the- 
ory introduces an extra friction term in the high-energy 
limit, which allows inflation to occur with much steeper 
potentials than usual [6J . This is favourable to get a 
second inflation today out of the same scalar field, i.e. 
a quintessence, from the behaviour that such potentials 
have at late times. 

At lower energies, the branewor ld-induced friction 
term is negligible and the scalar field spends a period 
completely dominated by its kinetic energy. In this era, 
usually called kination, the field density decreases very 
rapidly, p cx a~ 6 where a is the scale factor, permitting 
unusual reheating processes. The examples that will be 
explored later are reheating by gravitational particle pro- 
duction [7J, where radiation-like particles created during 
inflation can come to dominate the Universe during ki- 
nation, and curvaton reheating where energy density is 
stored in a second scalar field to be later released by its 
decay. 



The existence of a kination together with the possi- 
bility of having steep potentials makes this scenario ap- 
pealing for the implementation of quintessential inflation 
[H [5J. With a suitable choice of scalar field potential 
and reheating process one can try to fit observational 
data from inflation, big bang nucleosynthesis (BBN) and 
quintessence. There have been several studies of such 
models [5J IBHP2"]. but none of them satisfy all the con- 
straints brought by contemporary data. In this present 
work, we make a thorough review as well as expanding 
the space of models to fully understand to what extent 
it is possible to get inflation and quintessence out of the 
same scalar field in this simplified view of the braneworld 
paradigm. We will see that for the potentials analyzed, 
quintessential inflation cannot be achieved. 

In this paper, we discuss potentials with a simple 
structure, making only brief comments in Section IV on 
more complex ones, i.e. potentials that combine different 
forms at early and late times (for example Refs. [21 112)). 
Of course, we want potentials with good late-time be- 
haviours and, as will be shown, the exponential and in- 
verse power-law types are in that class. Therefore these 
are the two potentials analyzed by us in detail. 

The paper is organized as follows. In Section II we 
specify the models that we will test, presenting separately 
what is happening at early and late times. In Section III 
we impose constraints on these models based on their 
viability in comparison with observations. Again, the 
analyses of the constraints for the inflation period and 
late-time quintessence are treated independently. In Sec- 
tion IV we combine the constraints from both periods, to 
study the viability of quintessential inflation. Finally, in 
Section V we summarize the main results we obtained. 



II. BRANEWORLD QUINTESSENTIAL 
INFLATION MODELS 



The modified Fricdmann equation 
dimentional braneworld has the form [3] 
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where H is the Hubble parameter, A4 is a cosmological 
constant, is the four-dimensional Plank mass and A is 
the brane tension. The last term is related to the effect of 
bulk gravitons on the brane, but vanishes very quickly as 
inflation starts (spatial curvature is ignored for the same 
reason). Since we are interested in quintessence models, 
we assume that the four-dimensional cosmological con- 
stant is zero. 

The Friedmann equation then takes the simple form 
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where the second term will give an extra friction to the 
scalar-field evolution. The scalar field <j> that we will con- 
sider to give rise to inflation and quintessence is confined 
to the brane, following the standard evolution: 



3H<P 



dV 
' d(j>' 



(3) 



Inflation 



The simplest model of this class, presented in Ref. [6] , 
generates inflation with an exponential potential of the 
general form: 1 



V = Vq exp(— aip/rrii). 



(4) 



Here a is taken to be greater than vT&k, a regime in 
which the potential is too steep to sustain inflation in Ein- 
stein gravity. In the high-energy limit, with /> > A, the 
stronger friction term ensures that inflation can occur. 
During inflation, due to quantum fluctuations, radiation- 
like particles are created, in the generally-called gravita- 
tional particle production process [7J. At a low-energy 
limit, when p becomes smaller than A, the extra term of 
the Friedmann equation becomes unimportant and the 
field starts evolving as standard in a steep potential. In- 
flation ends and the field, experiencing such a steep po- 
tential, becomes dominated by its kinetic energy. Dur- 
ing this kination, the effect of the potential is negligi- 
ble and the energy density of the field decreases very 
rapidly as p oc a~ 6 . This abrupt loss of energy allows 
the gravitationally-produced radiation to dominate be- 
fore BBN. This reheats the Universe without requiring 
our scalar field to decay. 

The predictions of this model for the inflationary pa- 
rameters n s , the spectral index of scalar perturbations, 
and r, the ratio of tensor to scalar perturbations, can be 
easily obtained using the modified slow-roll parameters 



1 Note that many papers define a to be smaller by a factor 
by using the reduced Planck mass in the definition. 



of Ref. [5] , and have the simple form [B] 
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n s = 1 



N + l 
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where N is the number of e-foldings. It has been shown 
that in this case N can be unambiguously calculated to 
be 70 [§]. This results in the values n s — 0.944 and 
r = 0.33, which in combination are clearly excluded by 
current observations [T31 114) . 

This is one of the problems of this model. Another 
is related to a relic gravitational waves imprint. Grav- 
itational waves can be treated as a massless radiation 
field that decays as pow °c a_4 when it is in the sub- 
horizon regime. Short-scale gravitational waves, created 
by the same gravitational production mechanism, enter 
the subhorizon regime right after inflation, during the 
kinetic regime where the dominating background evolves 
as p oc a~ 6 . This means that these gravitational waves 
will have a boost in their relative energy density oc a 2 . 
Such a behaviour leaves an imprint that destroys BBN, 
excluding the model [9]. 

A solution to both these problems can come from the 
choice of another type of reheating, which brings the val- 
ues of \n s — 1| and r down as well as limiting the duration 
of the dominating kinetic period. This can be achieved 
with the inclusion of a curvaton field [HJ [16] whose de- 
cays reheat the Universe [TTJ [T7] . 

The curvaton is another field which coexists with the 
inflaton, that is subdominant and massless during infla- 
tion, and which can be responsible for the primordial 
curvature perturbations. For simplicity, we consider that 
this field has a quadratic potential of the form: 



U(a) = ^mV 



(7) 



where m is the mass of the curvaton. To ensure this is 
effectively massless during inflation we impose the condi- 
tion m <C Hf, the Hubble parameter at the end of infla- 
tion. Under these conditions, the value of the curvaton 
field, tr, remains constant during inflation |11[ 115]. apart 
from acquiring quantum-originated fluctuations from the 
expanding background. 

If we choose m correctly, during the kinetic period 
the field becomes effectively massive, i.e. m ~ H. This 
should happen at an early time to avoid the short-scale 
gravitational waves problem. The field then starts to os- 
cillate, which corresponds to behaviour as matter, and 
its energy density falls as p a oc a -3 [IS]. It then becomes 
dominant quite rapidly. 

The choice of m and of er,, the value of the field at 
the start of inflation, needs to be made in such a way 
that inflation and BBN can occur successfully. These 
constraints will be analyzed in detail in Section III. A. 

We will assume that the curvaton decays into radiation 
by some process with decay rate T, when T — H. The 
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field could start its decay while the kinetically-driven in- 
flaton is still dominating, which can bring some calcula- 
tional complications (TTJ [TS] . In this work, for simplicity, 
we will assume that the curvaton only starts decaying 
after being dominant. 

How are primordial perturbations created in this sce- 
nario? During inflation, the curvaton experiences quan- 
tum fluctuations, like the inflaton. These fluctuations of 
a freeze at horizon crossing and have a spectrum 
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where the star subscript denotes horizon crossing |15j . 
Then, after inflation, when the curvaton becomes mas- 
sive and oscillating, these fluctuations are converted into 
scalar curvature perturbations. The spectrum of these 
curvature perturbations has the form |15j 



J^Pl 



(9) 



The scale dependence of Pq is the same as that of P a 
and so the spectral index has the simple form 



n s = 1 - 2e, 



(10) 



where e is the slow-roll parameter. This matches the 
spectral index of primordial tensors and is usually differ- 
ent from that of the inflaton-induced perturbations. 

The parameter r (we follow the convention as in 
Ref. [H]), defined as 
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where /igw is the spectrum of tensor perturbations, also 
suffers alterations. Tensor perturbations in the metric are 
produced during inflation; they have the general form of 
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where F is a braneworld correction that in the high- 
energy limit takes the form 



F 
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In general, the scalar curvature perturbations Pq 
should be a sum of the perturbation of the inflaton and 
curvaton. In our case, we are assuming that the curvaton 
perturbations are dominant, so r has the form 



7271-0--= 



F 



(14) 



In Section III. A, where constraints are applied to the 
values of a and m, we will show how possible is the 
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FIG. 1: Inverse power-law potentials without curvaton. Com- 
parison of the values of n s and r with observations from 
WMAP five-year data, for different values of a. The expo- 
nential potential would lie under the a = 500 point. The plot 
shows the 68% and 95% confidence limits. The outer limits 
correspond to WMAP data alone, and the inner to WMAP 
data combined with supernovae and baryonic acoustic oscil- 
lation data. We can clearly see that this model is excluded 
by observations. Adapted from Ref. |13j . 



inclusion of a successful curvaton to the model, and if so 
how efficient it is in matching the WMAP observations. 

The other simple potential that we analyzed in our 
study is the inverse power law: 



V = V Q 



r/;,4 



(15) 



This potential allows an arbitrary amount of early- 
Universe inflation for any a > 2, as can be seen by 
the modified expression for the number of e-foldings in 
Ref. [5J. 

Just as in the case of the exponential potential, we 
could seek a good model by reheating the Universe with 
gravitational particle production. But the kinetic dom- 
inated period would still be too long, so the short-scale 
gravitational waves would remain a problem. 

Further, again using the modified slow-roll parameters 
IS], we can calculate n s and r: 



1 
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(16) 
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As can be seen in Fig. [T] these values are again excluded 
by WMAP. 
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So, as in the case of the exponential potential, to try to 
get a successful inflation, we need to include a curvaton 
in the model. In Section III. A we will show constraints on 
<Ji and m for this specific potential, to try to understand 
if a curvaton can in fact help the creation of inflation 
with good observational fit. 



B. Quintessence 

In the previous subsection, we sought to create a model 
that could give rise to successful inflation. Now it is nec- 
essary to check if these models can give rise to quint- 
essence, and if so, under which conditions. To do so, we 
need to analyze the late-time behaviours of the potentials 
studied. 

Kination won't last forever. The strong decrease of 
the energy density makes the kinetic energy decrease to 
a point where the potential is no longer negligible. The 
behaviours of such scalar fields that have been subdomi- 
nant for a certain period of time is well known [201 121] . 

Particularly interesting solutions when p$ < pback have 
the general name of scaling and tracking solutions, since 
in these cases the scalar field will 'follow' the evolution 
of the dominating background fluid. These two types of 
solution differ in the way the scalar field evolves. If we 
consider that the background fluid evolves as pback oc a~ n 
and that the scalar field evolves as p^ oc a _m , scaling so- 
lutions will be the ones where n = m and tracking the 
ones where n 7^ to. It is easy to see that such evolution 
is interesting, especially if we could get the scalar field 
tracking radiation and matter and only start dominating 
today. As shown in Ref. [5D], the only potentials that 
give rise to the above tracking behaviours are exponen- 
tials, inverse power laws and some positive power laws, 
although the last option is not interesting for connecting 
with inflation. Furthermore, it was also shown that these 
solutions are attractors. 

Scalar fields with exponential potentials have scaling 
solutions at late times, with n — m. Since the ratio 
between dominant fluid and scalar field is always a con- 
stant, there is no possibility for late-time quintessence if 
we are in the scaling regime. 

Scalar fields with inverse power laws have tracking 
solutions where n < to. This means that the scalar 
field energy density is increasing with time compared to 
the one of the dominant fluid. This is exactly the right 
condition to try to get quintessence. When the scalar 
field density approaches the one of the background 
fluid, the assumption that p^ < pback breaks down, and 
the field starts dominating, until its density is fl^, « 1 
and its equation of state is w ~ — 1. In Section III.B 
we will analyze whether, beside being able to produce 
acceleration at late times, this model can also satisfy 
the constraints on such acceleration today, i.e. ft fa 0.75 
and w < -0.8 P2QI]. 



nation period lasts for so long that the field overshoots 
its tracking solution and then freezes due to redshifting 
of the kinetic term, later approaching that solution from 
below. If the overshoot is so great that the field only ap- 
proaches the tracking solution around the epoch where 
tracking solution would be coming to dominate the fluid, 
we can enter directly into quintessence domination with- 
out a tracking period. This regime is referred to as creep- 
ing or thawing [52] quintessence, and was first explored 
in our present context in Ref. [8]. The energy density 
is initially approximately constant in this circumstance 
and the field behaves as a cosmological constant with 
w ps —1. When the energy density of the scalar field gets 
close enough to the background energy density, the field 
begins to evolve to w > — 1. 

With an inverse power-law, once thawing acceleration 
begins it will typically persevere. By contrast, the late- 
time solution for a steep exponential potential is a non- 
accelerating scaling solution, and in that case the field 
can only temporarily drive acceleration |23j . 

In Section III.B the above possibilities will be explored 
to see how well they can reproduce today's conditions on 
quintessence. 

III. CONSTRAINTS 
A. Inflation 

In this subsection we present the constraints that need 
to be imposed on the scalar field plus curvaton scenario, 
so that inflation occurs in a successful way. In fact, this 
will correspond to constraints on the space of values al- 
lowed for <7j and to, leaving as free parameters only N, 
the number of e-foldings, and a, the potential parameter. 

Before presenting the constraints, let us consider some 
useful relations. In the high-energy limit, 

H oc V (18) 

and V(j)/V( e ), where (i) and (e) refer to the beginning 
and end of inflation, can be expressed as a function of N, 
by using the modified expression for N given in Ref. [5]. 
For exponential potentials: 

-j^-=N + l (19) 
V {e) 

and for inverse power-law potentials: 

=N + l ^ ^ = (N + l)«/(«-2). (20) 

V(e } r {e) v [e) 

Using the modified expression for the slow-roll param- 
eter e, and setting it to unity, we can get an expression 
for V e . For exponential potentials we obtain 



An alternative regime occurs when the scalar field ki- 



(21) 
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and for inverse power-law potentials: 



V, 
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We don't so far know the value of <f>i e \ that shows up in 
this expression. To obtain it, we performed a numerical 
calculation of the equations of motion of the system dur- 
ing inflation. As a result, we obtained an approximate 
expression for Vu-\ depending only on V , a an d A: 



V (e) -0.06A Q /( Q - 2 )V 2/(Q " 2) ( 



(23) 



In the following calculations we will assume this expres- 
sion for Vr e ) . For simplicity, we also assume that kination 
starts immediately after inflation. 

Now we present the five conditions that must be 
satisfied for successful inflation: 

1. The curvaton needs to be subdominant when it 
becomes massive, to avoid a curvaton-driven inflation. 
This condition is expressed as: 



< 1 



(24) 



(m) 



where the subscript (m) refers to the period when the 
curvaton becomes massive. During kination, the Hubble 
parameter evolves as H oc a~ 3 , so 



m 



"kin 



and the Friedmann equation can be simplified to 
H 2 ^ 



3m j' 



We can then work out condition (|24 1 
Pi 



(m) 2p (a (k) /a (m) ) 



6m 2 



< 1 



(25) 



(26) 
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where we use (k) to refer to the beginning of kination. 
The condition is then simply: 



where the subscript '(eq)' refers to the curvaton-scalar 
field equivalence epoch. We can write 
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which is equivalent to having 
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So condition (30 1 is just 
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(31) 



(32) 
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4. Inflation perturbations cannot be important since 
we wish to get perturbations entirely from the curvaton. 
Without the presence of the curvaton, we could estimate 
the value of A no a from the COBE normalization of the 
density perturbations originated by inflation. To guaran- 
tee that curvaton-induced perturbations are much bigger 
than the ones originated by the inflaton, we just need to 
ensure that 



A < A r 



(34) 



5. The last condition that needs to be imposed is that 
short-scale gravitational waves are within limits. We can 
express this as 



Pgw 

Po 



Pgw 
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< 1. 



(35) 



(eq) 



Using the expression for the evolution of energy den- 
sity of gravitational waves during the kinetic period from 
Ref. [U, 



Pgw 
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(36) 



plus the amplitude of primordial gravitational waves, 
Eq. ( 12 ) which can be rewritten as 



af 3 

-V < IT- 
mi 4-7T 



(28) 



2. The curvaton needs to be subdominant at the end 
of inflation as well, to avoid a curvaton-driven inflation. 
This can be expressed as 



2V, 



< 1. 



(e) 



(29) 



3. The decay of the curvaton needs to occur before 
nucleosynthesis, so that this is successful, and after cur- 
vaton domination. We have the conditions: 



clco 



Kr 4 V 4 < r < h, 



(30) 
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(37) 



we can estimate our condition. 



Now, we want to evaluate these conditions for our two 
potentials. Conditions 1 and 3 are the same, as they 
don't depend on the shape of the potential. 

For exponential scalar field potentials, using expres- 
sions (19 1 and (21 1, condition 2 can be worked out as 
follows: 



m 2 cr 2 47r _ a 2 m 2 af(N + l) 2 



2a 2 A 



6H (i) m 



< 1. 



(38) 
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FIG. 2: Exponential potential with curvaton, showing constraints on the values of a and m for different a. From left to right, 
a — 5, a — 50 and a = 500. The viable space of parameters corresponds to the shaded area. The solid/orange line represents 
condition 1, the dotted/red line condition 2, the dashed/cyan line condition 3, the dash-dotted/green line condition 4 and the 
long-dashed/blue line condition 5. N was set to 70. The change in a doesn't strongly affect the results. 



So, we can write down the expression as 

547r 2 P c m! 



m 2 < 



a 2 (iV + l) 2 



(39) 



where P^ is the observationally-known spectrum of the 
Bardeen potential, which has the value of P^ = 2 x 10~ 9 
at the horizon scale. 

In condition 4, we need to insert the value that A would 
take if there was no curvaton in the model. This can be 
taken from the COBE normalization and looks like [5]: 



2.3 x 10- 10 (87i-) 3 m| 
a 6 {N + 1) 4 ' 



(40) 



A can be expressed as a combination of di and known 
parameters by using expressions ^ and ( 19 ) 
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Condition 4 is then just 



a? <C 



1 x 10 9 m| 
P c a 2 (iV+ l) 2 



(41) 
(42) 

(43) 



Condition 5, when expressions @, @, @ and 



(37) are plugged in, looks like: 



*(cg) 
a(fc) 



< 1 



(44) 



which, using Eqs. (21) and (31), can be written as (here 
we use the assumption that Hha = Hr e \) 

2/3 



/ 3rri4 



^a'(N+l)9P c af ^ 
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< 1. 
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FIG. 3: Exponential potential with curvaton. Comparison 
of the values of n a and r with observations from the WMAP 
five-year data. The data shown is as in Fig.[l] The parameter 
r can take different values depending on a. The interval of 
possible values of r is the same for all values of a. We can see 
that this model fits data successfully. Adapted from Ref . [13] . 



Substituting in the expression for A, Eq. (42), and rear- 
ranging, the condition is simply 



m > 



5 x 10 2 

r 3/2 



7714 



V(iv+i)P c V 



(46) 



Figure [2] shows these constraints in the m-a plane. 
There is a dependence on a, but as can be seen it is 
very weak. There is a comfortable parameter space 
where viable inflation can occur. In Fig. [3] we plotted 
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FIG. 4: Inverse power-law potential with curvaton, showing constraints on the values of a and m for different a. From left 
to right, a = 5, a = 50 and a = 500. The viable space of parameters corresponds to the shaded area. The solid/orange line 
represents condition 1, the dotted/red line condition 2, the dashed/cyan line condition 3, the dash-dotted/green line condition 
4 and the long-dashed/blue line condition 5. N was set to 70 and Vo to 10 _100 m4. We can see how the change in a, keeping 
Vo fixed, strongly affects the results. 



the result for n s and r for this case. The parameter n s is 
independent of a and a. The parameter r depends on a 
but the interval of possible values it can take is the same 
for each value of a. This model of inflation is in very 
good agreement with observations. In the next section 
we will see whether it can also produce a satisfactory 
quintessence. 

For inverse power-law potentials, we use always the 



approximate expression for Vu\, Eq. (23). Using this 
function, and following the steps as for hfq. (41), A can 
be expressed as 



A = \* Xl ^3 mla l V t7 2) \ . (47) 



16a 4 (iV+ l)2a/(a-2) 



So, condition 2 is 
m 2 < 0.12 



/ 2 x 10 3 7rm|P c 



\(N + l) 2 «/( a - 2 ) 



(48) 

For condition 4 we need to know what the value of A 
would be without the presence of the curvaton. It looks 
like [9] 



4a-8 



Pr 



(N + l) 4 



»-2 3a 

"+ 4 /47rA q+4 

' A . (49) 
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Using expressions (23) and (47), the condition can be 
written as 



ov < 6 x 10" 



!^5a 2 +10a-16^: 



2a-4 2a(a-2) p -2a+4 



(N + x^- 2a2 + 8tt + 16 7r 2a2 + 4 «+ 8 Q ,- 2 " 2 - 4a - 32 ^ 2 ( Q - 2 >(°+ 4 > 

(50) 



At this stage, Vq is a free parameter, which will be fixed 
when we impose that we want quintessence to start dom- 
ination today. However, Vo is not completely free, since 
the viability of the model depends strongly on the com- 
bination of this parameter and a. This can be seen by 
looking at our condition. Increasing Vo (with all other 
variables fixed) increases the right-hand side of the in- 
equality. Schematically, this corresponds to the dash-dot 
line of Fig. [4] moving towards higher values of a, increas- 
ing the space of parameters available. On the other hand, 
if both Vo and a were really small, the model would be 
unable to create good inflation. In Figs. [3] and [5j Vb is 
set to 10~ 100 m4 because it is a good initial guess for ac- 
celeration to happen today. This issue will be analyzed 
in more detail next section when we impose the correct 
value on Vo based on quintessence arguments. 

To get condition 5, we just follow the exact same steps 
as in Eqs. ([44]), |45| and ([46], leading to 
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(51) 



Figure [4] shows these five conditions for the inverse 
power-law potential in the m-cr space. In Fig. [5]we com- 
pare the range of values that n s and r can take with 
observations. We can see that the result is in good agree- 
ment with observations, although, for very small a, some 
values of a will give rise to too large an r. We can then 
impose a new constraint on the combination of the val- 
ues of a, Vo and a based on observations: r needs to be 
smaller than approximately 0.2 [T3") . 
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FIG. 5: Inverse power-law potential with curvaton. Com- 
parison of the values of n a and r with observations from the 
WMAP five-year data, for different a. The data shown is as 
in Fig. [I] For small a we can see that the parameter r can be 
too high compared to the observations. We can then impose 
an extra constraint on the combination of a, Vo and a. We 
can see that this model presents a big space of parameters 
that fits data successfully. Adapted from Ref. |13| . 



B. Quintessence 



In this section we will constrain the parameter space 
based on the observational data we have for quintessence. 
This is simply that today the dark energy density is ap- 
proximately fi^ f=a 0.75 and its equation of state has 
w < —0.8 [151 [13] . A successful model of quintessential 
inflation is then one for which the conditions at the end 
of inflation allow the scalar field to have these character- 
istics today. In other words, we want the scalar field to 
leave inflation at such a point in the potential and with 
such a velocity, that it will enter in a regime, being it 
tracking or thawing, that will resolve in a satisfactory 
domination today. 

Generally, smaller field values and larger velocities will 
increase the duration of kination and can force the scalar 
field to enter a frozen regime that can only redominatc 
through a thawing process. On the other hand, for the 
case of no velocity, the field always finds its tracking solu- 
tion, it being an attractor. For each a, the free variables 
that we can tune in order to change the conditions at 
the end of inflation are A and V). These completely de- 
termine <p e and <p e as we can see by setting the slow-roll 
parameter e to 1: in the standard formalism, e = 1 de- 
termines </> e [5] and, in the Hamilton- Jacobi formalism, 
€h = 1 determines cb e For example, for an inverse 




FIG. 6: and w for different a's for the tracking solution 
of an inverse power law. We can see that for a > 2, as it is 
required for the system to innate, it is impossible to match 
the observational constraints. 
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(52) 



(53) 



By choosing 4> e and </>, for each a, we can calculate 
the w of the scalar field today. We did this by perform- 
ing a numerical integration of the equations of motion of 
the system from the end of inflation until today. In our 
calculation we fixed = 0.75 today. Given a certain 
set of initial conditions, the program had the freedom to 
rescale Vq as many times as necessary to ensure that to- 
day indeed = 0.75. Once this condition was satisfied, 
one of the outputs of the calculation was w. 

The integration was done with respect to the number 
of e-foldings, so instead of Eq. (53) we use 



dN 



d<t> dt _ 2 
dt dN ~ a 



(54) 



Conveniently both (f>, 
ratio A/Vq and on a. 



and d(j)/dN depend only on the 



To check if any tracking solution would give a good 
quintessence behaviour, and since we already know that 
for the exponential case it doesn't, we run the calcula- 
tion for inverse power-law potentials with d<p/dN at end 
of inflation equal to zero. For any <p, because the solu- 
tion is an attractor, the scalar field found the tracking 
regime. The conclusion is that no tracking solution, for 
any a, could create a quintessence domination with the 
characteristic that we observe. The results are presented 
in Fig. [6] and match those of Ref. |S]. In particular, for 
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very big a's, the inverse power-law potential approaches 
the exponential one, and the scalar field never manages 
to dominate. Very small a like 1, which could give a vi- 
able quintessence, aren't interesting for us since we know 
that the system doesn't inflate for such a potential. 

To understand the true behaviour of the system, we 
run the calculation with the initial conditions as in the 



expressions (52) and (54). By choosing the ratio X/Vq, 



since the program will give the value of Vq , we are fixing 
A. We run it for inverse power laws, and we understand 
what occurs for exponential potentials by analyzing the 
large a limit. 

Firstly, we concluded that for A/Vb < 1 the model 
doesn't work because the scalar field never stops domi- 
nating. 

We also see that for small a's and small X/Vq the scalar 
field freezes and dominates at late times through a thaw- 
ing process. This only happens for a < 12. However, 
only for a < 6 we can actually achieve an equation of 
state w < —0.8 today. For greater a's, w is too big. 

As a increases it becomes harder for the field to have 
enough speed to freeze until domination by thawing, and 
it finds the tracking solution before dominating. Figure 
[7] shows two examples for the field thawing and finding 
its tracking solution before domination, respectively. 

For very large a, as expected, the field never manages 
to dominate at late times. We already knew that once 
the field finds its tracking solution, the model fails. 

We conclude, then, that for our conditions at the end 
of inflation, only for very particular cases of very small 
a and X/Vq are the observational constraints of quint- 
essence satisfied. Now we need to check if these con- 
straints are in agreement with the ones from inflation. 



IV. QUINTESSENTIAL INFLATION 

A satisfactory quintessence, as a manifestation of the 
inflaton at late times, is only achieved, as mentioned, for 
inverse power law potentials with 2 < a < 6 with low 
values of X/Vq (though bigger than 1). But for the full 
model to work, these constraints must be in agreement 
with inflation itself. 

When we fix X/Vq in the calculation, since Vq is an out- 
put value, we are fixing A. There are two tests we can do 
to check if inflation works successfully in the cases where 
quintessence does. First we can substitute the output 
Vq in the conditions for inflation and create a plot as in 
Figs. [2] and [4] to see if there is still a region of suitable 
parameters. If there is, the second test we can make is 



to substitute the value of A in expression (47) and see if 
the value of <Ji lays in such a region. 



In our case, Vq is of the order of 10 



for all the 



possible a's. This is a very low value. It implies that 
inflation would end at an energy scale of 10 _120 m|, which 
is much lower than the energy scale for nucleosynthesis. 
This incompatibility is clear in our plot of the parameter 
space. Figure [8] shows how by implementing the real 




FIG. 7: Examples of different evolutions of the scalar field. In 
the upper plot, a = 5. We see that the field is initially frozen 
and starts evolving until domination. This model is in good 
agreement with observations. In the lower plot, a = 20. We 
see that the field again begins in a frozen regime, but during 
its evolution, it finds the tracking solution. In this case, there 
is a disagreement with observations. In both plots, X/Vo = 1. 



value of Vq for a = 6, the parameter space disappears. 
This fact is even more drastic for smaller values of a. 

Hence we conclude that we cannot get quintessential 
inflation for the potentials studied. This is our main 
conclusion. 

Now, we can discuss how this result would change for 
other potentials. In this paper, we chose to study the 
simplest potentials that have interesting late-time be- 
haviours. More complex potentials with more freedom 
of parameters, when picked correctly, can be set to fit 
observations. In particular, if the potential is composed 
of two different regimes at early and late times, as in 
Ref. [5], it is always possible to fit the observations, just 
by choosing good potentials and by setting exactly where 
the late-time behaviour should start. Since our objective 
was to try to motivate the value of fi^ w 0.75 today by 
the conditions of the Universe at the end of inflation, 
such models aren't in the scope of our analysis. 

To end our discussion, although quintessential infla- 
tion is impossible to achieve in our framework, we want 
to highlight the fact that it is still possible to get in- 
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FIG. 8: Parameter space for a — 6 and Vo = W~ 120 m,p. 
We see that there is no possibility of successful inflation (the 
allowed parameter space is to the left of the dash-dotted/green 
line, above the dashed/cyan line and below the dotted/red 
line). For smaller values of a the discrepancy increases. 



flation. In such a case, for both potentials, the system 
would have to have a curvaton to give rise to the density 
perturbations and satisfy the constraints shown in Figs. [2] 
and [4] Additionally, to prevent it coming to redominatc 
too early, the inflaton would have to decay completely 
through some sort of reheating process, for instance a hy- 
brid inflation mechanism, leading to considerable further 
modelling complications. In this case, another explana- 
tion of the dark energy would need to be sought. 



obvious benefits in simplifying the structure of the cosmo- 
logical model, but must face much stronger observational 
challenges than models where separate components are 
responsible for each accelerated period. 

The range of models we have discussed is quite exten- 
sive. Having adopted the braneworld in order to have a 
natural end to inflation which does not require a feature 
in the potential, we then considered two substantially 
different reheating mechanisms — gravitational particle 
production and curvaton reheating — as well as two dif- 
ferent regimes of quintessence, tracking and thawing. In 
practice, gravitational particle production fails both to 
match the primordial spectra measured by WMAP and 
to give success nucleosynthesis due to excessive short- 
scale gravitational waves, forcing us to the curvaton case. 
Likewise, the conditions that inflation enforces on the 
potential immediately rule out the tracking models, as 
either the field fails to dominate at all at late times, or 
does so with an equation of state too far from —1. We 
are thus forced to the thawing regime, as first pointed 
out in Ref. [5]. It is then possible to find models capable 
of satisfying either the inflationary part of the picture, 
or the quintessence part. Sadly, however, we have found 
no models capable of satisfying both, the combined ob- 
servational constraints being failed by a large margin. 

Given the comprehensiveness of the models we have 
considered, we feel that current observations leave no 
room for quintessential inflation from simple poten- 
tials, with models assembled from the various known 
mechanisms of inflation, reheating, and quintessence. 
Quintessential inflation does, of course, remain possible 
provided one manipulates the potential into just the re- 
quired shape to bring the different regimes together, as in 
the original proposal of Peebles and Vilenkin [2] . To do so 
does not however appear any more attractive than main- 
taining two separate sectors, nor does it seem amenable 
to future observational tests. 



V. CONCLUSIONS 

We have presented a comprehensive series of models 
seeking to uncover the possibility of combining early Uni- 
verse inflation and present-day quintessence into a single 
field with a simple potential. Such a model would give 
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